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Abstract 

We study the delocalisation transition which takes places in one- 
dimensional disordered systems when the random potential exhibits 
specific long-range correlations. We consider the case of weak disor- 
der; using a systematic perturbative approach, we show how the de- 
localisation transition brings about a change of the scaling law of the 
inverse localisation length which ceases to be a quadratic function of 
the disorder strength and assumes a quartic form when the threshold 
separating the localised phase from the extended one is crossed. 

PACS numbers: 73.20. Jc, 72.15.Rn, 72.10.Bg 

Recently, one-dimensional (ID) models with correlated disorder have be- 
come the object of intense scrutiny, as a consequence of the discovery that 
specific long-range correlations of the random potential can create a contin- 
uum of extended states and give rise to mobility edges . The emergence 
of a phase of extended states implies that, when the disorder exhibits appro- 
priate long-range spatial correlations, even ID models can display a metal- 
insulator transition analogous to the Anderson transition which takes place in 
three-dimensional disordered systems. Although the results of the works [|I|, |2j 
were obtained for discrete lattices, it was soon found that long-range cor- 
relations can produce analogous phenomena also in continuous disordered 
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models ||. The mobility edge predicted in received an experimental con- 
firmation in where the phenomenon was observed in the transmission of 
microwaves in a single-mode guide with a random array of correlated scat- 
terers. Further progress along this line of research would make possible the 
realization of ID devices with specified mobility edges which could be used 
as window filters in electronic, acoustic or photonic structures, as suggested 
in §. 

The theoretical and practical importance of the delocalisation transition 
in ID models makes desirable to obtain a complete understanding of the role 
played by correlations in the formation of mobility edges. Unfortunately, 
the theoretical comprehension of the phenomenon achieved in and in the 
successive analytical works is not totally satisfactory, because it rests on the 
assumption of weak disorder and on the use of a second-order approxima- 
tion scheme. Thus, the delocalisation transition has been defined by the 
condition that the second-order term of the inverse localisation length (or 
Lyapunov exponent) vanish over a continuum interval of electronic energies, 
while higher-order terms of the Lyapunov exponent have been completely 
neglected. This simplifying approach is effective but incomplete, because it 
does not clarify the true nature of the "extended" states, leaving open the 
question of whether they are completely delocalised or not. The aim of this 
work is to shed light on this point and thus provide a better understanding 
of the delocalisation transition in ID models. To this end, we make use of 
a systematic perturbative approach that allows us to compute the fourth- 
order term of the Lyapunov exponent. Armed with this result, we are able 
to ascertain the nature of the extended states and to discuss the differences 
between Gaussian and non-Gaussian potentials which appear in the localisa- 
tion properties of ID models at this refined level of description. 

In what follows, we consider the Hamiltonian continuous model defined 
by the Schrodinger equation 

- (h 2 /2m)i)"(x) + eU(x)ip(x) = Eip{x) (1) 

which describes an electron of energy E moving under the influence of a 
random potential U(x). We study the case of weak disorder, identified by 
the condition e — > 0. To simplify the form of the equations, in the rest of this 
paper we adopt a system of units such that h 2 /2m = 1. We consider electrons 
of positive energy, so that we can write E = k 2 ; since the wavevector k enters 
model (|1]) only via the energy, we can restrict our attention to the case k > 
without loss of generality. 
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We define the statistical properties of the disorder through the moments 
of the potential U(x). We assume that the model under study is spatially 
homogeneous in the mean; as a consequence, the n-th moment of the poten- 
tial 

Xn(xi, x 2 ,..., x n -i) = (U (x)U(x + Xx) . . .U(x + X n -i)) 

depends only on n — 1 relative coordinates. Here and in the following we use 
angular brackets (. . .) to denote the average over disorder realizations. We 
also want the average features of model ([[]) to be invariant under a change 
of sign of the potential; therefore we assume that all moments of odd order 
vanish. Finally, for the sake of mathematical convenience, we temporarily 
suppose that a correlation length l c can be defined such that the values U(x\) 
and U(x2) of the potential are statistically independent when the distance 
between the points x\ and x 2 is much larger than l c . To study the case of 
long-range correlations, we will subsequently relax this hypothesis by taking 
the limit l c — > oo. 

To determine the average behaviour of the eigenstates of Eq. (|T|), we 
introduce the inverse localisation length 

X(k) = Jim i- ln(^ 2 (x)k 2 + ^' 2 (x)). (2) 
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As a first step to compute the Lyapunov exponent (0), we derive a differential 
equation for ip 2 (x) and ip' 2 (x). Using Eq. (|l|) as a starting point, it is easy 
to show that the vector 

ip (x) 
u{x) = I i)' 2 {x)/k 2 
^(x)ip'(x)/k 



{A + et{x)B)u, (3) 



obeys the equation 

du 
dx 

where ^(x) = U (x) / k 2 is the scaled potential and the symbols A and B stand 
for the matrices 



A 



/ 2k 

-2k I and B 

V -k k 



/ 

V -k 




At this point our problem consists in extracting information on the average 
vector (u(x)} from the stochastic equation for u(x); this can be done using 
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a method developed by Van Kampen ||, which replaces Eq. (|3|) with an 
ordinary differential equation of the form 

d(u(x)) , .. . . . 

dx () 

The generator K(x) in Eq. ([|) is a sure operator that can be expressed in 
terms of a cumulant expansion 

oo 

K(x) = J> n K n (x) (5) 

n=0 

where the partial generators K n (x) are functions of specific combinations 
(known as "ordered cumulants" ) of the moments of the random matrix 

M(x) = £(x) exp (-Ax) Bexp (Ax) . (6) 

Van Kampen has established a set of systematic rules for constructing all 
the terms of the series (|5|); applying his prescriptions one obtains that all 
generators of odd order are zero because the odd moments of the potential 
vanish. As for the even-order generators, the zero-th order term has the 
simple form K = A, which corresponds to ignoring the random potential. 
The first term where disorder manifests its effect is the second-order generator 

K 2 (x) = / dxi exp (Ax) (M(x)M(xi)) exp (-Ax) , 
Jo 

which is proportional to the two-point correlator; a refinement of this result 
is obtained by taking into account the fourth-order term 

K 4 (x) = Jq dxi Jq 1 dx2 Jq 2 dx3 exp (—Ax) 

x [(M(x)M(x 1 )M(x 2 )M(x 3 )) 

-(M(x)M(x 1 ))(M(x 2 )M(x 3 )) 

-(M(x)M(x 2 ))(M(x 1 )M(x 3 )) 

- (M(x)M(x 3 ))(M(xi)M(x 2 ))] exp (Ax) 

whose integrand constitutes a non-trivial example of ordered cumulant of the 
matrix (^j). Although the rules of Van Kampen make possible to derive every 
partial generator K n (x), for the purpose of this work we can truncate the 
series (|5|) and neglect all terms of higher order than the fourth. As can be 
seen from Eq. (|j), the behaviour of (u(x)) for x — > oo is determined by the 
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eigenvalue with the largest real part of the asymptotic generator K(oo) = 
lim^oo K(x); we are thus led to the conclusion that, in order to determine 
the Lyapunov exponent (0) with fourth-order accuracy, we have to solve the 
secular equation for the truncated asymptotic generator 

K(oo) = A + e 2 K 2 (oo) +e 4 K 4 (oo). (7) 

Carrying out the calculations, one obtains that the inverse localisation 
length (Q) can be written as 

X(k) = e 2 X 2 (k) + e A (xf\k) + X[ NG \k)) + o (e 4 ) (8) 

where the second-order term has the standard form 

1 f°° 

^2{k) = —— X2(x)cos(2kx)dx. (9) 

4fiT Jo 

As for the fourth-order term of the Lyapunov exponent (§|), its first compo- 
nent can be written as 



A 4 G) (^) = W¥>»(M) bc(M) +4^(0,0)] 

-gp [<f c (k, 0) + 2y2 c (0, 0)] J °° X2(x)x cos(2kx)dx 

+s^ ( Ps(k, 0) X2(x)x sm(2kx)dx 

+IF Jo°° bc(0, x)ip c (k, x) - ip 2 c {0, x) cos(2A;a;)] dx 

where the functions <f c {k, x) and <f s {k, x) are defined as 

<Pc(k,x) = J™ X2(y) cos(2ky)dy 
<Ps(k,x) = 4°° x%{v) sin(2ky)dy ' 

The second component turns out to be 

Xf G) (k) = pr / °° dx dy J™ dzA 4 (y, x, z) 
x [cos(2%) cos{2kz — 2kx) — cos(2kz)} 



(10) 



(11) 



where the symbol A 4 (xi, x 2 , £3) represents the fourth cumulant of the random 
potential. These formulae show that, in contrast to the relative simplicity 
of the second-order result (|9|), the fourth-order terms ( |i~0"|) and (|TT|) exhibit 
a non-trivial dependence on the binary correlator and on the fourth cumu- 
lant, respectively. The greater complexity of the results is a consequence 
of the more sophisticated description of localisation, which is visualised as 
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an interference effect generated by double scatterings of the electron in the 
second-order scheme, whereas the fourth-order approximation also takes into 
account quadruple scattering processes. The more accurate description of the 
localisation processes makes possible to draw a distinction between Gaussian 
and non-Gaussian disorders. These two classes of random potentials are nec- 
essarily undifferentiated within the framework of the second-order approxi- 
mation, since the second-order Lyapunov exponent @ depends only on the 
second moment of the potential; in the fourth-order description, however, the 
differences between the two kinds of potentials emerge in the component (fTT[) 
which, being a function of the cumulant, vanishes if, and only if, the disorder 
is Gaussian. In conclusion, the separation of the fourth-order term of the 
Lyapunov exponent in the two components ( PD and (|TT|) is justified by the 
distinct physical character of these two parts: while the former is common 
to all Gaussian and non-Gaussian potentials with the same two-point cor- 
relation function, the latter describes the specific effect of the non-Gaussian 
nature of the disorder. 

To analyse the case of disorder with long-range correlations, we evaluated 
the form taken by expressions @ and (|10D when the correlation function is 

Xa(s) = a 2 Sin(2M exp(-x/Q (12) 
x 

and we subsequently took the limit l c — > oo, which reduces the correlation 
function to the form 

sin(2fc c x) 

X2{X) = 0-2 • 13 

x 

This approach is generally successful, but it fails for k = k c because for this 
value of the wavevector the asymptotic limit of the generator (|5|) exists only 
as long as l c is finite. Our results for the Lyapunov exponent, therefore, are 
not valid in a small neighbourhood of k c . 

Inserting the correlation function flT2|) in Eq. (|9|) and taking the limit 
l c — > oo, one obtains that the second-order Lyapunov exponent has the form 

^> ={?[:: l<v K • (i4 » 

Eq. (|14|) shows that long-range correlations of the form (|T3|) make the second- 
order Lyapunov exponent vanish when the wavevector k exceeds a critical 
value k c . This result was interpreted in |2[] as evidence that long-range cor- 
relation of the potential can create a continuum of extended states, with a 
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mobility edge at k = k c . Second-order theory, however, cannot precise the 
nature of the "extended" states nor define their spatial extension. To over- 
come this limit, one can compute the fourth-order term ([II]) of the Lyapunov 
exponent. Evaluating expression (|T0|) for the correlation function ([Tj|) and 
then taking the limit l c oo, one arrives at the result 



\f\k) 




< k < k c 
2 k 



where the functions Ai(fc) and A 2 (k) are defined as 



A 2 (k) = 



A. 

1A: 4 



5?r i / fc+fc e V 
32fc 111 \k-k c ) 

8k 111 \k-k c ) 



7T fc 2 +fcfc e + fcg 
8 k c (kl~k' 2 ) 

T 2k c —k 



I(k) 
8 k c (k-k c ) I\k) 



with 



/(*) 



si (2k c x) cos(2kx)dx. 



(15) 



The symbol si (a;) in the previous equation represents the sine integral defined 
as si(x) = — dtsm(t)/t. In Fig. [I] we show the plot of the fourth-order 
Lyapunov exponent ([15]) as a function of the wavevector k. As already ob- 



Figure 1: Fourth-order term X^ik)/^ vs. k/k, 
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served, the divergence for k = k c must be disregarded because the asymptotic 
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generator (|7]) is not denned for this value of the wavevector in the limit of 
long-range correlations. In the interval < k < k c the second-order Lya- 
punov exponent ( |H|) is positive and the fourth-order term (|i5|) constitutes 
just a small correction which enhances or decreases the localisation of elec- 
tronic states according to whether it assumes positive or negative values. 
The main interest of Eq. (p~5|) , however, lies in the fact that it shows that 
the fourth-order Lyapunov exponent is positive in the interval k c < k < 2k c , 



where the second-order term (|14]) vanishes. This result clarifies the nature of 
the mobility edge at k = k c : the electronic states are exponentially localised 
on both sides of the threshold, but the localisation length is increased by a 
factor 0(l/e 2 ) when the wavevector exceeds the critical value k c . For weak 
disorder, the corresponding change in the electron mobility is so conspicu- 
ous that one can legitimately speak of "delocalisation transition", even if 
the "extended" states are not completely delocalised like Bloch states in a 



crystal lattice. A second noteworthy aspect of the fourth-order term ( |15| ) 
is that it vanishes for k > 2k c so that, in the case of Gaussian disorder, a 
second transition takes place at k — 2k c with the electrons becoming even 
more delocalised. 

So far our discussion has been centred on the term (|H3), which allows a 
complete analysis of Gaussian potentials. For non-Gaussian disorder, how- 
ever, one must also take into account the term (|TT|) which enhances localisa- 
tion. As an example, one can consider the case of a non-Gaussian potential 
with long-range binary correlations of the form (|T3j) and with an exponen- 
tially decaying fourth cumulant 

A 4 (x,y,z) = (7 4 exp[-/3(|x| + \y\ + \z\)] . (16) 

In this case the fourth-order term of the Lyapunov exponent is the sum of 
the component (|I5|) and of the additional term 



4 1 ; 8k 2 {(3 2 + k 2 ){(3 2 + 4k 2 ){9f3 2 + 4k 2 ) 1 ; 

which is obtained by inserting expression (|16D in Eq. (0). Since the extra 
term ( |TTD is always positive, all electronic states with k > k c are exponentially 
localised with inverse localisation length proportional to the fourth power of 
the disorder strength, X(k) oc e 4 . 

The main conclusion of this study is that the delocalisation transition 
occurring in ID models with long-range correlated disorder consists in a 
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change of the scaling behaviour of the Lyapunov exponent when the electronic 
energy crosses a critical threshold: more specifically, the Lyapunov exponent, 
which is a quadratic function of the disorder strength in the localised phase, 
assumes a quartic form in the extended phase. Another relevant remark is 
that non-Gaussian potentials produce a stronger localisation effect than their 
Gaussian counterparts. 

The author wishes to thank S. Ruffo and A. Politi for useful discussions 
and is grateful to F. M. Izrailev for valuable comments and suggestions. 
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